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^~^An  approximate  analytical  model  has  been  developed  for 
computing  the  first  cycle  of  the  shock  structure  in  plumes  from 
underexpanded  axisymmetric  supersonic  nozzles.  The  physical  basis 
of  the  model  is  the  same  as  that  of  the  Abbett  model,  i.e.,  the 
Mach  disk  location  is  determined  by  the  coupling  of  the  flow  in 
the  vicinity  of  the  Mach  disk  with  the  reaccelerating  transonic 
flow  downstream  of  the  disk.  -However,  instead  ef  the  Abbett — 
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-technique  of  — 1-teratlve  flow-field  ooapu  tat  lens  between  the  trial 
Jlaoh  disk-  and  the  transonie  region,  £he  present  model  is  based 
on  a local  compatibility  condition  for  the  flow  just  downstream 
of  the  Mach  disk.  To  evaluate  the  local  compatibility  condition, 
the  generalized  shape  of  the  slip  line  which  originates  at  the 
shock  triple  point  must  be  known.  The  Abbett  method  was  used  to 
compute  slip  line  geometries  for  a number  of  plume  flows; 
correlation  equations  for  the  slip  line  geometrical  parameters 
were  then  developed..  The  flow  field  upstream  of  the  Mach  disk 
is  computed  with  the  rotational  method  of  character is tioe . At 
eaett  point  ©n-the  boundary  shock  wave,  a trial  triple  point  is 
assumed -and  the- -loeal— compatibility  condition— is  evaluated;  the 
Maoh  dioJels  located  wtrere  the  compatibility  equation  is  satisfied. 
After  the  Mach  disk  is  located,  the  supersonic  flow  outside  of 
the  slip  line  is  computed  for  the  prescribed  slip  line  geometry. 
The  approximate  analytical  model  requires  relatively  small 
computation  times  and  yields  acceptably  accurate  flow-field 
solutions  over  the  range  of  conditions  for  which  the  slip  line 
correlation  equations  are  valid. 
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PREFACE 

The  work  reported  herein  was  conducted  by  the  Arnold  Engineering 
Development  Center  (AEDC),  Air  Force  Systems  Command  (AFSC),  under 
Program  Element  65807F.  The  results  were  obtained  by  AR0,  Inc.  (a 
subsidiary  of  Sverdrup  & Parcel  and  Associates,  Inc.),  contract  operator 
of  AEDC,  AFSC,  Arnold  Air  Force  Station,  Tennessee.  The  work  was  done 
under  ARO  Project  Nos.  RF423,  R33P-60A,  and  R33A-02A.  The  authors  of 
this  report  were  C.  E.  Peters  and  W.  J.  Phares,  ARO,  Inc. 

The  Fox-Abbett  computations  which  are  discussed  in  this  report 
were  carried  out  by  J.  H.  Fox  of  the  Central  Data  Processing  Division 
of  ARO,  Inc.  The  manuscript  (ARO  Control  No.  ARO-ETF-TR- 76-85)  was 
submitted  for  publication  on  August  5,  1976. 
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1.0  INTRODUCTION 

Numerous  experimental  investigations  (e.g.,  Refs.  1 through  4)  have 
defined  the  general  structure  of  the  axisymnetric  plume  which  is  formed 
when  the  flow  from  an  underexpanded  sonic  or  supersonic  nozzle  exhausts 
into  a quiescent  ambient  medium  (Fig.  1).  A boundary  shock  wave  is 
formed  by  the  coalescence  of  compression  waves  from  tlie  curved  jet  bound- 
ary.  As  the  boundary  shock  wave  propagates  downstream,  it  becomes 
stronger  and  curves  toward  the  axis.  For  slightly  underexpanded  nozzle 
operation  (pn/pa>  < 1.5),  the  boundary  shock  wave  appears  to  undergo  a 
regular  reflection  at  the  axis.  As  Pn/poD  is  increased,  however,  the 
boundary  shock  no  longer  undergoes  a regular  reflection  at  the  axis; 
instead,  a nearly  normal  shock  wave,  called  the  "Mach  disk"  or  "Riemann 
wave",  occurs  (Fig.  1).  The  boundary  shock,  the  Mach  disk  and  the  re- 
flected shock  are  joined  at  a triple  point. 


Figure  1.  Multi  pie  shock  wave  structure  in  plumes. 

For  moderately  underexpanded  nozzle  operation  (pn/pa>  up  to  10  or 

20),  the  plume  shock  structure  has  been  observed  to  repeat  with  increasing 

distance  (Fig.  1).  Thus,  the  flow  behind  the  Mach  disk  reaccelerates  to 

supersonic  speeds  then  passes  through  another  shock  system.  The  observed 

number  of  repetitions  of  the  shock  system  depends  on  p /p  , but  the 

n 
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inviscid  shock  wave  structure  is  usually  obliterated  after  two  or  three 
cycles  by  the  turbulent  mixing  layer  between  the  Jet  and  the  ambient  fluid. 

If  Pn/po  is  very  large,  perhaps  100,  then  the  maximum  diameter  of 
the  plume  and  the  diameter  of  the  Mach  disk  are  much  larger  than  the 
nozzle  exit  diameter.  It  is  not  clear  from  the  available  experiments 
whether  the  plume  shock  structure  is  repetitive  at  such  high  pressure 
ratios.  Almost  all  of  the  experimental  information  on  plume  shock 
structure  has  been  obtained  by  optical  methods;  for  highly  underexpanded 
plumes,  the  shock  structure  downstream  of  the  first  Mach  disk  is  obscured 
by  the  turbulent  mixing  along  the  plume  boundary. 


Various  approximate  techniques  have  been  proposed  for  predicting 
the  location  of  the  first  Mach  disk.  The  empirical  correlation  equation 
of  Lewis  and  Carlson  (Ref.  5)  yields  acceptably  accurate  Mach  disk 
positions,  but  the  technique  provides  no  information  about  the  plume 
structure  other  than  the  Mach  disk  location.  Adamson  and  Nicholls 
(Ref.  2)  assumed  that  the  Mach  disk  occurs  at  a position  where  the 
pressure  behind  the  normal  shock  on  the  centerline  is  exactly  equal 
to  the  ambient  pressure,  p^.  This  method  requires  only  that  the  flow 
properties  along  the  plume  centerline  be  defined.  The  Mach  disk  loca- 
tions predicted  with  the  Adamson  and  Nicholls  method  are,  in  general, 
less  accurate  than  those  predicted  with  the  Lewis  and  Carlson  empirical 
equation.  No  information  is  provided  about  the  two-dimensional  distri- 
butions of  flow  properties  in  the  plume  with  either  method. 

Several  investigators  have  proposed  models  for  Mach  disk  location 
which  are  based  on  detailed  computation  of  the  plume  flow  field.  The 
inviscid  flow  upstream  of  the  Mach  disk  can  be  computed  with  the  rota- 
tional method  of  characteristics.  However,  the  method  of  characteristics 
offers  no  insight  into  the  mechanism  of  Mach  disk  formation;  the  flow 
field,  including  the  boundary  shock  wave,  can  be  computed  downstream 
beyond  the  experimental  location  of  the  Mach  disk,  with  no  indication 


6 


AEOC-TR -76-129 


that  the  solution  really  does  not  occur.  Therefore,  in  addition  to  the 
rotational  method  of  characteristics,  the  "field"  methods  require  assump- 
tions about  the  mechanism  of  Mach  disk  formation. 

Eastman  and  Radtke  (Ref.  6)  assumed  that  the  Mach  disk  occurs  at 
a position  where  the  static  pressure  behind  the  computed  boundary  shock 
reaches  a minimum.  Bowyer,  et  al.,  (Ref.  7)  assumed  that  the  Mach  disk 
is  a normal  shock  at  the  triple  point.  Therefore,  the  Mach  disk  is 
predicted  to  occur  at  the  station  where  the  usual  triple  point  conditions 
(P3  ■ p^,  63  ■ 0^)  are  satisfied  and  where  (Fig.  2). 


The  Eastman  and  Radtke  and  Bowyer,  et  al.,  field  methods  are  attrac- 
tive because  they  provide  a prediction  of  the  entire  plume  flow  field 
upstream  of  the  Mach  disk.  But  these  methods  predict  Mach  disk  locations 
which  are  in  only  fair  agreement  with  experiment  (Ref.  8).  Furthermore, 
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the  assumed  mechanisms  of  Mach  disk  formation  do  not  have  a sound  physical 
basis. 

The  "field"  methods  were  generally  unsuccessful  until  a physically 
perceptive  mechanism  for  Mach  disk  formation  was  postulated.  The 
correct  physical  mechanism  was  pointed  out  by  Peters  (Ref.  9),  Peters 
and  Phares  (Ref.  8),  and  Abbett  (Ref.  10):  the  Mach  disk  forms  such 

that  its  position  is  compatible  with  the  reaccelerating  transonic  flow 
downstream  of  the  disk.  Abbett  was  the  first  to  develop  a fairly 
rigorous  analysis  of  Mach  disk  formation.  The  flow  upstream  of  the 
Mach  disk  is  computed  with  the  rotational  method  of  characteristics. 
Downstream  of  the  Mach  disk,  the  flow  outside  of  the  slip  line  (Fig.  1) 
is  also  computed  with  the  method  of  characteristics.  The  slip  line  is 
assumed  to  be  inviscid  and  the  flow  between  the  slip  line  and  the  axis 
is  assumed  to  be  one-dimensional  and  isentropic.  A trial  position  of 
the  Mach  disk  is  assumed,  and  the  combined  supersonic-subsonic  flow 
downstream  of  the  disk  is  computed.  If  the  Mach  disk  is  assumed  too 
far  upstream,  then  the  slip  line  will  reach  a minimum  radius,  a "throat", 
before  the  flow  on  the  axis  reaches  sonic  speed.  If  the  Mach  disk  is 
assumed  too  far  downstream,  then  the  flow  on  the  axis  will  become  sonic 
before  the  slip  line  radius  reaches  a minimum.  The  position  of  the 
Mach  disk  is  iterated  until  the  flow  on  the  axis  becomes  sonic  at  a 
station  where  the  slip  line  forms  a throat. 

The  throat-like  condition  downstream  of  the  Mach  disk  is  a saddle 
point  singularity  in  the  flow-field  solution.  As  usual,  with  saddle- 
point  singularities,  it  is  not  too  difficult  to  bracket  the  position  of 
the  Mach  disk  to  acceptable  uncertainty.  But  it  is  very  difficult  to 
establish  the  Mach  disk  position  precisely  enough  so  that  the  downstream 
solution  will  pass  smoothly  through  the  transonic  region.  Even  if 
the  Mach  disk  position  is  bracketed  very  tightly,  to  one  part  in  10^ 
or  10^,  the  solution  will  in  general  not  pass  smoothly  through  the 
transonic  region.  However,  it  is  possible  to  use  extrapolation  techniques 
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to  approximate  the  solution  through  the  transonic  region.  Once  the 
plume  flow  field  downstream  of  the  Mach  disk  is  fully  supersonic,  the 
Abbett  method  can,  in  principle,  be  used  to  compute  the  second  cycle  of 
the  plume  shock  structure.  However,  the  assumption  of  one-dimensional 
flow  behind  the  Mach  disk  may  be  less  valid  for  the  second  Mach  disk 
than  for  the  first. 

Abbett  presented  the  solution  for  only  one  plume  flow  field.  However, 
Fox,  who  developed  a plume  shock  structure  analysis  which  is  essentially 
the  same  as  Abbett* s,  has  demonstrated  that  the  method  yields  accurate 
predictions  of  the  Mach  disk  location  for  a variety  of  nozzle  exit  Mach 
numbers  and  pressure  ratios  between  3 and  16  (Ref.  11). 

Using  the  concept  that  the  Mach  disk  is  uniquely  related  to  the 
transonic  flow  downstream  of  the  disk,  Peters  and  Phares  (Ref.  8)  devel- 
oped a plume  shock  structure  analysis  which  is  much  simpler  than  Abbett* s. 
Instead  of  iteratively  calculating  the  flow  field  downstream  to  the  throat 
region,  the  generalized  shape  of  the  slip  line  downstream  of  the  triple 
point  was  assumed.  A third  triple  point  boundary  condition  (in  addition 
to  p^  - p^  and  0^  ••  0^)  was  written: 

dp/dx|3  - dp/dx|^  (1) 

That  is,  the  slip  line  just  downstream  of  the  triple  point  must  allow 
the  flow  on  either  side  of  the  slip  line  to  develop  compatibly.  The 
flow  upstream  of  the  Mach  disk  is  computed  with  the  rotational  method 
of  characteristics.  At  each  axial  location  in  the  boundary  shock  compu- 
tation, a triple  point  computation  is  made  to  define  p^  and  0^  (Fig.  2). 

By  assuming  that  the  slip  line  is  a circular  arc,  dp/dx|^  can  be  computed, 

based  on  only  the  local  properties.  Linearized  supersonic  flow  theory 
is  used  to  compute  dp/dx|3<  The  Mach  disk  position  is  established  when 
Eq.  (1)  is  satisfied. 
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Although  based  on  a sound  physical  mechanism,  the  Peters  and 
Phares  model  yields  only  fair  predictions  of  the  experimental  Mach  disk 
position  (Ref.  8).  The  failure  of  the  model  can  be  attributed  to  two 
causes.  First,  the  flow  outside  of  the  slip  line  (region  3 of  Fig.  2) 
is  highly  rotational  and  the  assumed  adequacy  of  linearized  theory  is 
incorrect.  Second,  the  assumed  circular  arc  slip  line,  from  the  triple 
point  to  the  sonic  station,  is  incorrect.  Experiments  (e.g.,  Ref.  12) 
have  shown  that  the  slip  line  radius  downstream  of  the  triple  point  first 
increases  to  a maximum  before  turning  toward  the  axis.  The  Abbett 
analysis  predicts  such  a complex  slip  line  geometry. 

The  analysis  described  in  this  report  is  basically  similar  to  that 
described  in  Ref.  8,  in  that  the  coupling  between  the  Mach  disk  and 
the  transonic  region  is  described  in  terms  of  the  third  triple  point 
boundary  condition  (Eq.  1).  However,  the  description  of  the  slip 
line  geometry  has  been  improved,  and  the  assumption  of  linearized 
supersonic  flow  in  region  3 has  been  discarded.  Physically  realistic 
results  are  obtained  without  the  time-consuming  flow-field  iterations  of 
the  Abbett  analysis. 

2.0  DEVELOPMENT  OF  ANALYTICAL  MODEL 

The  approach  taken  in  the  present  analysis  has  been  to  use  the 
Fox  version  of  the  Abbett  analysis  to  compute  the  slip  line  geometries 
for  a series  of  plume  flows.  "Empirical"  equations  were  then  developed 
to  correlate  the  results  of  the  Abbett  analysis  for  the  slip  line  geo- 
metrical parameters.  The  empirical  slip  line  parameters  are  used,  in 
conjunction  with  a single  rotational  method  of  characteristics  field 
point  calculation  in  region  3 of  Fig.  2,  to  evaluate  the  validity  of 
each  trial  position  of  the  triple  point. 


The  slip  line  geometries  for  18  plume  flows  were  computed  with 
the  Fox-Abbett  computer  program.  The  18  flows  represented  a random 
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sampling  for  M^  from  1.0  to  2.5,  Pn/P00  from  3 to  16,  y from  1.2  to  1.4, 

and  nozzle  exit  half-angles  from  0 to  15  deg.  (In  all  cases,  the  nozzle 

exit  flow  was  assumed  to  be  source-like).  Two  characteristic  length 

scales  for  the  slip  line  were  identified  (Fig.  3):  (1)  the  distance, 

from  the  Mach  disk  to  the  station  where  the  slip  line  radius  is 

a maximum,  r and  (2)  the  distance,  from  the  Mach  disk  to  the  sonic 

smax  2 

point  on  the  axis. 


Figure  3.  Ratio  of  length  scales  for  the  slip  line. 


The  ratio  of  to  was  f°und  t0  correlate  well  as  a function  of 

the  Mach  disk  radius,  r /r  . As  shown  in  Fig.  3,  the  Fox-Abbett 

m n 

results  are  well  described  by 


H./i,  - 0.264  (r  /r  ) 
1 / m n 


-0.37 


(2) 
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The  ratio  yrffl  was  found  t0  correlate  well  with  the  flow  angle 
upstream  of  the  triple  point,  0r  Aa  shown  in  Fig.  4,  the  Fox-Abbett 
results  are  well  described  by 

yrn  - 130(ei)-1*55  (3) 


Equations  (2)  and  (3)  are  combined  to  yield 

1,/r  - 34. 3(r  /r  )”0’ 37(01)”1 ' 55 
I tn  is  n 1 


(4) 


The  slip  line  geometry  between  and  (xa  + y is  approximately 
parabolic  (Fig.  and  is  described  by 


(5) 
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But  tan  0.-2  (r  - r )/£.,  so  that 
4 smax  m 1 

1 2 

r * r + (tan  8.)(x  - x ) - -r  (tan  6.)(x  - x ) /l.  (c\ 

am  4 m 2 4 ml 

and 

r - r + 4 i.  tan  6.  n\ 

smax  m 2 1 4 \>  > 


Figure  5.  Parabolic  slip  line  geometry  between  x - x„,  and  x - xm  +8,. 


Downstream  of  the  maximum  slip  line  radius,  the  slip  line  shape  (Fig.  6) 
is  approximately  described  by 


(' 


x - x 


max 


x*  - x 


max 


■)] 


0.67 


(8) 
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Figure  6.  Slip  line  geometry  between  the  maximum  radius  station  and  the  sonic  station. 

Triple  point  computation  procedure  - The  flow  field  upstream  of  the 
Mach  disk  is  computed  with  the  rotational  method  of  characteristics. 

The  computer  program  is  an  extension  of  the  "2-D  Core  Theory"  of  Refs. 

13  and  14,  which  includes  the  capability  for  a fully  coupled  computation 
of  the  turbulent  mixing  layer  along  the  plume  boundary.  However,  in 
order  to  consistently  compare  the  present  model  with  the  inviscid  Fox- 
Abbett  model,  the  turbulent  mixing  at  the  plume  boundary  was  suppressed 
in  most  of  the  computations  presented  in  this  report. 

A trial  triple  point  is  assumed  at  each  axial  station  in  the  com- 
putation of  the  boundary  shock  wave.  The  Rankine-Hugoniot  equations 
are  used,  along  with  the  triple  point  boundary  conditions  (p^  “ p^, 

“ 0^),  to  compute  p^  and  0^.  The  values  of  and  0^  upstream  of 
the  trial  triple  point  are  used  to  compute  from  Eq.  (4).  Then  Eq.  (6) 
is  used  to  establish  the  geometry  of  the  slip  line.  With  the  known  slip 
line  geometry  and  the  assumption  that  the  subsonic  flow  is  one-dimensional 
and  lsentropic,  the  pressure  distribution  and  flow  angle  distribution 
along  the  slip  line  are  established. 
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From  a point  "a"  just  downstream  of  the  triple  point  (Fig.  7),  at 
which  the  pressure  and  flow  angle  are  known,  a left  running  characteristic 
is  generated  to  establish  the  coordinates  and  flow  properties  of  point  "b" 
just  behind  the  reflected  shock  wave.  Then  the  location  of  point  "c"  on 
the  slip  line  is  computed  with  the  equation  which  describes  the  geometry 
of  the  right  running  characteristic  between  points  "b"  and  "c";  the 
axial  distance  between  x and  point  "c"  is  typically  only  a few  percent 

IS 

of  . All  flow  properties  at  point  "c"  are  fully  defined  by  the 
equations  which  describe  the  flow  in  the  one- dimensional  streamtube. 

If  the  compatibility  equation  for  the  right  running  characteristic 
between  points  "b"  and  "c"  is  not  satisfied,  then  the  assumed  slip  line 
geometry  does  not  permit  consistent  development  of  the  flow  on  either 
side  of  the  slip  line,  and  Eq.  (1)  is  not  satisfied.  Therefore,  the 
computation  of  the  flow  field  is  continued  downstream,  with  a trial 
triple  point  assumed  at  the  next  axial  station.  At  some  axial  station, 
the  compatibility  equation  for  the  right  running  characteristic  will  be 
satisfied;  this  station  is  the  Mach  disk  location. 
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Computation  procedure  downstream  of  the  Mach  disk  - After  the  Mach 

disk  position  is  established,  the  flow-field  computation  can  be  continued 

downstream  beyond  the  sonic  station  on  the  axis.  Equations  (3)  and  (4) 

are  used  to  compute  i,  and  r is  computed  from  Eq.  (7),  and  r * 
r 1 2*  smax  s 

is  computed  from  the  isentropic  area  ratio  equation.  Equation  (6) 
describes  the  slip  line  geometry  between  x^  and  (x^  + and  Eq.  (8) 

describes  the  slip  line  geometry  between  (x^  + and  (x^  + 

With  the  inviscid  slip  line  geometry  prescribed,  the  rotational 
method  of  characteristics  is  used  to  compute  the  supersonic  flow  outside 
of  the  slip  line.  The  computation  is  continued  downstream  until  the 
static  pressure  on  the  slip  line  is  slightly  less  than  the  sonic  pressure 
for  the  one-dimensional  streamtube.  Because  the  supersonic  flow  field 
is  generated  along  a prescribed  geometry  (the  slip  line),  no  singularity 
appears  in  the  solution. 

3.0  DISCUSSION  OF  RESULTS 

The  Mach  disk  locations  predicted  with  the  approximate  plume  shock 
structure  model  are  compared  with  the  experimental  results  of  Love,  et  al. , 
(Ref.  3)  in  Fig.  8.  The  nozzle  air  flow  exhausted  into  quiescent 
surroundings.  All  of  the  results  to  be  discussed  were  obtained  with 
nozzles  which  were  designed  to  provide  uniform  and  parallel  flow  at 
the  exit  plane. 

Computed  and  experimental  Mach  disk  positions  are  shown  in  Fig.  8a 

for  Mn  « 2.0.  The  Mach  disk  positions  predicted  by  the  approximate 

analysis  are  about  ten-percent  upstream  of  the  experimental  positions 

at  the  lower  values  of  p /p  , but  only  a few  percent  upstream  of 

n • 

experiment  for  Pn/PB  “ 15.  The  Fox-Abbett  analysis  provides  excellent 
predictions  of  the  experimental  Mach  disk  locations. 

Computed  and  experimental  Mach  disk  locations  are  shown  in  Figs. 

8b,  c,  and  d for  Mr  of  1.0,  1.5,  and  2.5,  respectively.  In  addition  to 
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the  experimental  results  of  Love,  et  al.,  for  M - 2.5,  two  experimental 

n 

points  obtained  at  the  AEDC  for  * 2.44  are  shown  in  Fig.  8d.  In 

general,  the  predictions  of  the  approximate  model  are  no  more  than  ten 

percent  in  error  at  the  low-pressure  ratios,  and  the  predictions  improve 

as  p /p  increases  toward  16. 
n • 


1 2 4 6 8 10  12  14  16 


Pressure  Ratio, 
a.  M„  - 2.0 


Pressure  Ritlo,  pn^a 

b.  M„  - 1.0 

Figure  8.  Prediction  of  the  axial  location  of  Mach  disk. 
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Pressure  Ratio.  pn/pa 

c.  Mn  - 1.B 


Prassura  Ratio,  Pn^PoD 

d.  Mn  - 2.5 
Figure  8.  Concluded. 


Very  little  information  is  available  on  the  flow  structure  down 
stream  of  the  Mach  disk.  R.  A.  Paulk  and  C.  E.  Peters  conducted  a 
simple  experiment  to  define  the  sonic  station  downstream  of  the  Mach 
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disk  in  the  plume  from  a Mach  number  2.44  nozzle  operated  at  p /p  ■ 

tl  00 

10*7  (Ref.  11).  Blunt  bodies  of  various  sizes  were  placed  on  the  plume 
axis  downstream  of  the  Mach  disk.  The  bodies  were  then  translated 
upstream  toward  the  Mach  disk.  At  some  critical  distance  between  the 
Mach  disk  and  the  body,  the  Mach  disk  was  observed,  with  a shadowgraph, 
to  be  affected  by  the  body.  By  extrapolating  the  critical  distance 
function  to  zero  body  size  (a  negligible  disturbance),  the  axial  distance, 
<>2*  between  the  Mach  disk  and  the  sonic  station  was  established.  The 
experimental  results,  as  well  as  the  predictions  of  the  Fox-Abbett  and 
approximate  analyses,  are  shown  below. 


x /r 
m n 

Vrn 

Experiment 

13.1 

6.4 

Fox-Abbett  Analysis 

13.06 

6.44 

Approximate  Analysis 

12.35 

6.45 

Both  predicted  values  are  in  excellent  agreement  with  the  experimental 
value. 

In  the  computations  which  have  been  discussed  up  to  this  point,  the 
turbulent  mixing  at  the  plume  boundary  has  been  neglected.  Inclusion  of 
the  turbulent  mixing  causes  the  boundary  shock  wave  to  be  stronger  than 
in  the  inviscid  boundary  case.  As  a result,  the  Mach  disk  is  predicted 
to  occur  upstream  of  the  location  predicted  for  an  inviscid  plume  bound- 
ary. However,  as  shown  in  Fig.  8c  for  Mr  - 1.5,  inclusion  of  the  bound- 
ary mixing  changes  the  predicted  value  of  x by  only  about  five  percent. 

m 

The  results  which  are  presented  in  Figs.  8a  through  d indicate  that 
the  approximate  plume  model  yields  satisfactory  predictions  of  the  Mach 
disk  location  for  p / p between  3 and  16,  that  is,  satisfactory  results 
are  obtained  over  the  p^p^  range  which  was  used  to  develop  the  empiri- 
cal slip  line  geometrical  parameters.  However,  only  fair  results 
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are  usually  obtained  for  Pn/Pa>  values  of  1.5  to  2.5.  The  empirical 
slip  line  functions  must  be  refined  if  the  approximate  model  is  to 
yield  good  results  for  low-pressure  ratios.  The  validity  of  the  empiri- 
cal slip  line  functions  at  pressure  ratios  greater  than  16  was  not 
investigated. 

Because  the  iterative  flow-field  computations  of  the  Abbett  method 
have  been  eliminated  in  the  approximate  shock  structure  model,  rela- 
tively small  computation  times  are  required.  Typically,  computation 
of  the  first  cycle  of  a plume  requires  one  minute  on  an  IBM  370/165 
digital  computer.  The  computation  time  is  not  significantly  affected 
by  including  the  turbulent  mixing  along  the  plume  boundary. 

4.0  CONCLUDING  REMARKS 

The  main  objective  of  this  study,  to  develop  a computationally 
efficient  and  reasonably  accurate  analytical  model  for  the  first  cycle 
of  the  plume  shock  structure,  has  been  achieved.  However,  acceptably 
accurate  predictions  of  the  Mach  disk  location  can  be  expected  only 
for  Pn/P„,  between  3 and  16,  the  range  for  which  the  slip  line  geomet- 
rical correlations  were  developed.  These  correlations  must  be  improved 
if  the  Mach  disk  location  is  to  be  accurately  predicted  for  other 
pressure  ratios.  Until  improved  correlation  equations  are  available, 
it  is  suggested  that,  at  pressure  ratios  less  than  3,  the  Mach  disk 
location  be  established  with  the  Lewis  and  Carlson  empirical  equation. 
The  techniques  of  the  present  study  can  then  be  used  to  compute  the 
flow  through  the  transonic  region  downstream  of  the  Mach  disk. 

In  principle,  the  present  analytical  model,  as  well  as  the  Abbett 
model,  1b  applicable  to  (he  computation  of  the  second  cycle  of  the 
plume  shock  structure.  However,  the  flow  approaching  the  second  Mach 
disk  is  much  more  rotational  than  the  flow  approaching  the  first  disk. 
Therefore,  the  assumption  of  one-dlmenslonal  flow  behind  the  second 
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Mach  disk  may  not  be  acceptably  accurate.  In  addition,  computation 
of  the  second  cycle  should  probably  include  the  mixing  along  the  slip 
line  which  extends  downstream  from  the  first  Mach  disk.  Clearly,  the 
applicability  of  presently  available  models  to  the  computation  of  the 
second  cycle  of  the  plume  requires  further  study.  However,  further 
study  will  be  hampered  by  the  lack  of  detailed  experimental  information 
on  the  flow  structure  downstream  of  the  first  Mach  disk. 
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NOMENCLATURE 

^ Axial  distance  from  the  Mach  disk  to  the  station  where  the 

slip  line  radius  is  a maximum 

Axial  distance  from  the  Mach  disk  to  the  station  where  flow 
on  the  axis  is  sonic 

M Mach  number 

M Nozzle  exit  Mach  number 

n 

p Static  pressure 

p Nozzle  exit  static  pressure 

rn 

p^  Ambient  pressure 

r Mach  disk  radius 

m 

r Nozzle  exit  radius 

n 

r Slip  line  radius 

s 

r Maximum  slip  line  radius 

smax 

r * Slip  line  radius  at  the  sonic  station 

s 

x Axial  coordinate 

x Axial  distance  from  the  nozzle  exit  plane  to  the  Mach  disk 

m 


1 ' *8 
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x Axial  distance  from  the  nozzle  exit  plane  to  the  station 

max 

where  the  slip  line  radius  is  a maximum,  * x + it. 


x*  Axial  distance  from  the  nozzle  exit  plane  to  the  station 

where  the  flow  on  the  axis  is  sonic,  = x + 

m 2 


Y Ratio  of  specific  heats 


0 Flow  angle 


SUBSCRIPTS 


1,2, 3, 4 Flow  regions  defined  in  Fig.  2 
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